The w-mode instability of ultracompact relativistic stars 
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In this Letter we show that the, mainly spacetime, w-modes may become unstable in a rotating 
ultracompact (R < 3M) relativistic star. We provide results for the axial modes of a rotating star 
in the slow rotation approximation, and demonstrate that an instability may become active in stars 
with an ergosphere. These results suggest that the so-called "ergoregion instability" discussed by 
other authors is, in fact, associated with the w-modes of the star. The unstable modes are found 
to grow on a timescale of seconds to minutes. We discuss briefly whether the instability is likely to 
have any astrophysical relevance. 

(N : 

• Introduction. — Just like a rotating black hole, a spinning neutron star can develop an ergosphere. The ergosphere 
is defined as a region where all observers must be dragged along with the rotation, and no observer can remain at rest. 
£^ ■ Mathematically, the ergosphere is distinguished as the part of a stationary spacetime where a spacelike Killing vector 
becomes timelike. Hence, the boundary of the ergosphere is defined by the change in the sign of the [tt] component 
^vq ■ of the metric. In the case of relativistic stars, an ergosphere may develop inside in the stellar fluid if the star is 
\ sufficiently compact (R < 3M). 

It is by now well known that a neutron star spacetime with an ergosphere may become unstable. This was first 
demonstrated by Friedman ffl] for scalar and electromagnetic perturbations. The problem was studied in detail by 
Comins and Schutz Q who considered scalar fields in the background of a slowly rotating star. Their approach had the 
advantage that the problem could be reduced to one dimension, which facilitated the use of the WKB approximation 
' to calculate the unstable modes of oscillation (for large m, where the perturbations are assumed to be proportional 
' to e lm ^). Comins and Schutz found that the growth time of the instability would be very large (of the order of the 
age of the universe). This would mean that, even if sufficiently compact and rapidly rotating stars would form, the 
instability would not be astrophysically important. More recently, the problem was revisited by Yoshida and Eriguchi 
■ H who developed a new numerical scheme which allowed them to extended the results of Comins and Schutz to 
smaller values of m. Their results showed that the low-order (in m) modes would grow much faster than indicated 
Q ■ by the asymptotic results of Comins and Schutz. 

Since the first studies of the ergoregion instability were carried out the existence of pulsation modes that are 
essentially due to the spacetime itself has been established Q . These are known as the w-modes j| , and they only 
. . | rely on the background curvature generated by the fluid in order to exist. The actual coupling to the stellar fluid is 
weak, and the mode mainly excites oscillations of the spacetime. An extreme case concerns axial perturbations of a 
non-rotating perfect fluid star. In this case no non-trivial fluid modes of oscillations exist. Yet they have an infinite 
set of w-modes. A similar set of spacetime modes exists for polar perturbations, although in that case there are 
also non-trivial fluid modes (for example, the f-mode and the acoustic p-modes), see |g| for an exhaustive discussion. 
In general, the w-modes have high frequencies, with the fundamental mode in the range 5-12kHz, and very short 
damping times, of the order of a tenth of a millisecond. Before the w-modes were discovered Chandrashekhar and 
Ferrari j(| suggested that ultracompact stars (with radius smaller than 3M) can support oscillation modes due to the 
trapping of gravitational waves by the curvature potential (familiar from studies of perturbed black holes) . Detailed 
studies || Q have shown that these long-lived modes are, in fact, the first few w-modes of an ultracompact star. The 
key difference is that, because these modes are trapped inside the peak of the curvature potential their damping can 
be very slow (the damping time is roughly proportional to the width of the potential barrier) . 

As we will demonstrate below, these trapped w-modes may become unstable if the star is set into rotation. We 
will show that this instability only becomes active after the star has developed an ergosphere. This means that the 
gravitational-wave version of the ergoregion instability § 1 sets in through the spacetime w-modes. 

The w-mode instability. — In order to establish the existence of the w-mode instability we study the effect of rotation 
on the axial w-modes within the slow rotation formalism. In writing down the perturbation equations we will keep 
only terms of first order in the angular frequency of the star (Q) . Hence, we account for the relativistic frame dragging, 
but do not consider the rotational deformation of the star (which is an order f2 2 effect). In addition, we neglect the 
coupling between axial and polar perturbations. 
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If we focus on purely axial perturbations, the perturbed metric can be written as 

ds 2 = ds 2 + 2 (h l m (t, r)dt + h[ m (t, r)dr) (- sin" 1 dd^Y^dO + sin 6d e Y h 



(1) 



where Yi m = Yi m (8,(j>) denote the scalar spherical harmonics. The unperturbed metric ds^ represents a non-rotating 
stellar model with a first order rotational correction u(r) (the frame dragging) in the [i, 0]-component: 
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The fluid is assumed to rotate with uniform angular velocity f2, and the axial component of the fluid velocity pertur- 
bation can be expanded as 



(3) 
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With these definitions Einstein's field equations reduce to three equations for the three variables h m , h l ™ and U lm 
H . An equivalent set of equations in the ADM- formalism is given in |], [l(J . Assuming a harmonic time dependence 
e lat , the perturbation equations can be written as two ODEs for ho and hi and an algebraic relation for U ; 



/in 



A' 



ma/ (A - 2) 



A(er — muj) 
too/ 



hi - 


'2m 2 r 2 I 


A ( 





LU ,2 e- 2X 



1&ttw 2 (j> + e) 



A(fT — mu) 2mtu + A(a — mO) 



+ [a — muj) 



ih 



1 



r A(er — mto) 



h + 



e 2v {K - 2) 



= 4(^-A) 



U = Aire- 2 " {p + e) 



I - 



r 2 (a — mto) 
2mw 



muj ; 



ih% 



2mvo + A(cr — mQ) 



(4) 
(5) 
(6) 



where p is the pressure, e the energy density and A = 1(1 + 1). 

Equations ffl) and (^) become singular for certain values of a. This singularity is due to the existence of a continuous 
spectrum J^, [ll], [T^]. The existence of this continuous spectrum has no impact on our study of the w- modes, while it 
greatly affects the study of r-modes. 

Given the above equations it is worth noticing that in equation (^) the coefficient of ho includes terms which are 
formally of second order in the rotation rate tt. Within the slow rotation approximation these terms can be neglected. 
However, this would not be consistent for inertial modes (like the r-mode) since their frequencies are such that a ~ fi. 
Only in this way can the correct growth times of the r-modes be calculated ||. For the w- modes we have typically 
a ^> f2 and we can neglect second order terms in Q. However, we nevertheless use the above set of equations in order 
to be able to compute both r- and w-modes with the same numerical code. 




FIG. 1: Oscillation frequency (left frame) and damping/growth time (right frame) of the first w-mode of ultracompact stars. 
The rotation rate of the star is represented by e = £1/Qk- The w-mode instability is active beyond the point where a neutral 
mode exists, i.e. whenever the mode frequency is negative. As is clear from the data, the instability can grow on a timescale 
of a few tenths of a second to minutes in very rapidly rotating stars. 

Our numerical results for ultracompact stars are illustrated in Fig. |l|. The two frames show the oscillation frequency 
and damping/growth time of the first w-mode for a star with mass I.4M Q , the canonical value for a neutron star and 
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four different compactness ratios. The illustrated mode is such that its pattern moves backwards with respect to the 
rotation as fl — > 0. As the rotation rate, represented by e = Q,/Qk (where VLk is the mass shedding limit), increases 
the frequency of the mode decreases and passes through zero at a critical value. Beyond this point the emerging 
gravitational waves drive the mode unstable. The onset of instability is also signalled by a singularity in the growth 
time since a marginally unstable mode would not radiate at all. These singularities are obvious in the right frame of 
Fig. 0. 

It is relevant to compare the critical rotation rate at which the first w-mode becomes unstable to the Kepler limit. 
For the cases shown in Fig. |l| we find that the instability becomes active at 0.190^ for R/M = 2.26, at 0.33ilx 
for R/M = 2.28, at 0AAfl K for R/M = 2.30 and at 0.83fi K for R/M = 2.40. In other words, for stars near the 
absolute limit of compactness allowed by General Relativity the w-modes become unstable already at quite low rates 
of rotation (below 20% of the mass-shedding limit). As is clear for the right frame in Fig. [I] the associated growth 
times can be as short as a few tenths of seconds to minutes. Thus the growth time of the w-mode instability can be 
very short (Our results indicate that the unstable gravitational perturbations grow at least four orders of magnitude 
faster than the scalar field perturbations considered by eg Yoshida and Eriguchi S. Nevertheless, we find that the 
unstable w-modes grow slower than the I = m = 2 r-mode in all cases we have considered. We have also verified that 
the growth times scale with the rotation rate as ~ f2 2Z+2 as one would expect. 

The class of "trapped" w-modes which becomes unstable consists of a finite number of modes. The number of the 
"trapped" modes increases with the compactness of the star and several of these modes become unstable for higher 
rotational rates. For example for a star with compactness R/M = 2.26 near the mass shedding limit more than 4 
modes become unstable. 
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FIG. 2: The size of the ergosphere for our four stellar models is shown for different values of e. The ergosphere is present above 
a certain critical value of e. We also show (as straight lines) the critical value of e at which the first w-mode becomes unstable 
for each stellar model. As is clear from this data, the w-mode instability only becomes active after an ergosphere has formed. 

We now want to establish that the w-modes will only be unstable in stars that have developed an ergosphere. To 
test the existence of an ergosphere we follow Schutz and Comins p3| and approximate gtt as 

g tt = ^ + lo 2 t 2 sin 2 6 (7) 

Given this, an ergosphere is present whenever gtt > 0, i.e. restricting ourselves to the equatorial plane 9 = n/2 we 
require 

e"<™. (8) 

Figure § illustrates the ergosphere (in the equatorial plane) for our ultracompact stellar models and varying rates 
of rotation. Also shown in the figure are the critical values of the rotation parameter e at which the first w-mode 
becomes unstable. As is clear from the figure, the w-mode only becomes unstable once the star has developed an 
ergosphere. 

Discussion — In this Letter we have demonstrated that the spacetime w-modes may become unstable in ultracompact 
rotating stars. Our results show that the instability come into operation after the star has developed an ergosphere, 
and that the unstable modes can grow quite fast — on a timescale of seconds to minutes. The w-mode instability 
is a variation of the gravitational-wave driven (CFS) instability, which was discovered by Chandrasekhar jl4j for the 
Maclaurin spheroids and subsequently investigated in detail by Friedman and Schutz [|l5] . The key difference between 
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the w-mode instability and the standard CFS instability of fluid oscillation modes (like the axial r-mode |Q ) is that 
the w-modes are mainly due to spacetime oscillations. They are gravitational waves trapped (in the present case) 
inside the peak of the curvature potential. Their instability can be understood from the fact that an oscillation mode 
can be associated with a negative energy provided that the star has an ergosphere. 

As an interesting aside it is worth commenting on the fact that despite them being endowed with an ergosphere, 
rotating black holes do not suffer an analogous instability. The reason for this is simple: In the black hole case, the 
leakage of energy through the horizon proceeds fast enough to stabilize the system. 

The demonstration that the w-modes may become unstable is interesting from a conceptual point of view. However, 
the astrophysical relevance of this instability is questionable. The main reason for this is that it is very difficult to 
construct realistic stable stellar models which would exhibit an ergosphere (Tt|, Indeed, there are very few 

"realistic" supranuclear equations of state that permit stable stellar models significantly more compact than R w 3M. 
A possibility that cannot yet be completely ruled out is that of ultracompact strange stars. One can, in principle, 
generate very compact quark stars but we have not yet any observational evidence for the existence of such objects. 

Finally, it is worth discussing two features that make the w-mode instability differ somewhat from the more familiar 
CFS instability of fluid oscillations modes. The first concerns dissipation. Since the w-modes are spacetime modes 
that hardly excite any fluid motion there are no apparent dissipation mechanisms. The viscous damping of the 
instability will certainly be significantly weaker than that of (say) the r-modes since the viscosity will only act on a 
relatively small part of the w-mode energy. The second point concerns the saturation of the unstable mode. Recent 
work suggests that the unstable r-modes saturate at a relatively low amplitude because of strong coupling to short 
wavelength inertial modes. The weak coupling to the stellar fluid in the present case may mean that this mechanism 
will be significantly less efficient, if of any relevance, in the case of unstable w-modes. In fact, one could argue that 
these modes ought to saturate because of nonlinear spacetime effects. More detailed investigations into these issues 
would be of interest since the may help improve our understanding of nonlinear dynamical spacetimes. 
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